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Abstract 

We determine the form factor expansion of the one-point functions in integrable 
quantum field theory at finite temperature and find that it is simpler than previously 
conjectured. We show that no singularities are left in the final expression provided 
that the operator is local with respect to the particles and argue that the divergences 
arising in the non-local case are related to the absence of spontaneous symmetry 
breaking on the cylinder. As a specific application, we give the first terms of the low 
temperature expansion of the one-point functions for the Ising model in a magnetic 
field. 
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Quantum field theory at finite temperature is a subject of both theoretical and experi- 
mental interest. In the euclidean formulation, it corresponds to work with an immaginary 
time compactified on a circle whose circumference R coincides with the inverse tempera- 
ture. The average of a generic product X of local operators is given by 

(x) - TlXe ~ RH m 

{ }R ~ Tre- RH ' { } 

where H is the Hamiltonian. 

In two dimensions, a non-perturbative study of the properties of the finite tempera- 
ture case should be achievable exploiting the results of integrable quantum field theory. In 
particular, a natural idea is that of approaching the computation of the correlation func- 
tions in a way similar to that successfully used in the zero temperature case, and evaluate 
([!]) on the basis of the multi-particle asymptotic states of the zero temperature theory. 
At the end, the problem should be reduced to summing over the matrix elements of the 
operators between the vacuum and the n-particle asymptotic states (the form factors) 
which are known exactly in integrable theories. 

The one-point functions have a non-trivial temperature dependence and provide the 
first test that this form factor method needs to pass. Following the analysis of the free 
fermionic case |I[ , an extension to the generic integrable case exploiting also the results of 
the thermodynamic Bethe ansatz (TBA) was conjectured in Ref. 0. More precisely, it was 
proposed that the 'pseudoenergy' given by the TBA should substitute the ordinary energy 
of the particles in the thermal weight function entering the form factor sum. Lukyanov 
found agreement with this conjecture in his semi-classical study of the sinh-Gordon model 
at finite temperature ||. 

In this note, we analyse explicitely the form factor expansion of the one-point functions 
in finite temperature integrable theories and find that the final expression, Eq. (|I6|) , is 
simpler than proposed in the past: the weight function is the same than in the free cases 
and all the effects of the interaction are contained in the form factors. Along the way we 
discuss a number of delicate issues originating from the fact that the expression (|I|) is a 
sum over matrix elements between identical asymptotic states which contain singularities. 
We explain why and in which cases the final result is free of divergences. 

This analysis also allows a clarification of the status of the conjecture of Ref. 0. In 
fact, the expression for the one-point functions proposed in that work differs from our 
Eq. (|16D not only in the thermal weight, but also in the operator dependent part. In 
Refs. H [|, the latter does not coincide with the form factors (|l^) but is such that its 
combination with the modified weight function reproduces the TBA result for the one- 
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point function of the trace of the stress-energy tensor. At present it is not clear how a 
similar expression involving the pseudoenergies could be derived for other operators. In 
any case, due to the complicated temperature dependence of the pseudoenergies, it would 
not be the explicit low temperature expansion (|16|) . 

Consider an integrable quantum field theory whose spectrum contains a single neutral 
particle of mass M. As usual, let us use the rapidity variable 9 to parameterise the energy 
and momentum of a particle as (p°, p 1 ) = (M cosh 9, M sinh 9). Then relativistic invariant 
quantities depend on the rapidity differences only. We denote S(9i — 9 2 ) the scattering 
amplitude of two particles with rapidities 9i and 9i- It satisfies the crossing symmetry 
relation 

S[9) = S{m - 9) . (2) 

Integrability shows up in two respects f|. First, it forbids any inelastic process, so that 
the unitarity condition takes the simple form 

S(9)S(-9) = 1 . (3) 

Second, it induces the complete factorisation of multi-particle scattering amplitudes into 
the product of the two-particle ones. 

Take now a scalar operator in this theory which is local with respect to the 

particles, and denote its matrix elements on the asymptotic multi-particle states as 

i£ B (C . . . , e[\e u . . . , 9 n ) = (C • • • , 0[mo)\Oi, • • • A) • (4) 

The contribution of the n-particle asymptotic state to Tr e~ HR is 

fn(R) = ~ } T^w [dOx... d9 n F* n (9 n , . . . , 9^9,, . . . , 9 n ) e~ E " R , (5) 
n! (2ir) n J 

where 

n 

E n = Mj2cosh9i , (6) 

i=l 

is the total energy of the asymptotic state, and we used the property ^(x) = e lP ^ x ^^(0)e~ lPtJ 
f®{R) is of order e~ nMR for MR — > oo. After defining 

oo 

/*(*) = E£(*)> (7) 

n=0 

we can write 

<*>« = £g , (8) 
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with / denoting the identity operator. 

It is our goal to rewrite this rather implicit expression containing the matrix elements 
F® n as a sum over the form factors F® n which are exactly computable in integrable field 
theories || Technically, this problem has several points of contact with Ref. |§, where 
non-integrable field theories were considered as perturbations of the integrable ones. In 
that context the matrix elements F* n determine the first order correction to the scattering 
process m —>■ n under the action of the perturbing operator <£>. 

The reduction to the form factors F * n is performed iterating m times the crossing 
relation g, [| 

F m,n(®'mi • • • > ^5.1^1) ■ ■ ■ , 9 n ) = F^_ ln+1 (9' m) . . . , 9 2 \9[ + 17T,9 1 , . . . , 9 n ) + 
n i—1 

2tt £ 5(9[ - Bi) II S(9 k - 9[) F*_ l n _ 1 (^, . . . , 9' 2 \9 Xl . . . , . . . , fl B ) . (9) 

i=i fe=i 

The second term in the r.h.s. accounts for the disconnected parts that appear if the crossed 
particle hits a particle with exactly the same momentum and annihilates it. It contains 
the product of the scattering amplitudes with the particles that need to be crossed on the 
way. We recall in fact that the relation 

I • • • 9i, 9 i+ \ . . .) — S(9i — 9 i+ i) \ . . . 9 i+ i, Bi. . .) , (10) 

is characteristic of integrable field theories. The situation of collinearity responsible for 
the disconnected parts also induces an 'annihilation' pole in the matrix elements with 
residue Rl 



ires ei= g 2+i7T F (\9i, . . . , 9 n ) 



F* (|0 S ,...A). (11) 



l-\{S{9 l -9 2 ) 

i=i 

Since (||) contains matrix elements involving two identical sets of momenta, it is obvious 
that we will have to care about both annihilation poles and singular disconnected parts 
containing 5(0) factors. We will see in a moment, however, that these potential sources 
of trouble drop out and leave a final result which is perfectly finite. 
Using @ the first few /*(-R) are found to be 

/* = ^ = (0|$|0) = ($) fl=oo , 

fi(R) = ^ / d9, e^ R [J* + 2nS{0)jf\ , 

f!(R) = \ J d9,d9 2 e~ E * R [F*{9 Xt 9 2 ) + *k{5{0 x - 9 2 )S(0) + 5(0)) 

+ (2nnS(0)5 2 (9 1 - 9 2 ) + 5 2 (0))^], 
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f!(R) = ^ j^y 3 J d9 l( i9 2 d6 3 e~ E * R [J^(B 1 , 9 2 , 9 3 ) + 6^23(0)5(9, - 9 3 ) + 5(0))Ff(9 u 6 2 ) 

+ 3(2n) 2 (25(9 1 - 92)5(0! - 9 3 ) + S(0)5 2 (9 1 - 9 2 ) + 23(0)5(0)5(9, - 9 2 ) + 5 2 (0))F* 
+ (2tt) 3 (2(J(0 1 - 9 2 )5(9, - 9 3 )5(9 2 - 9 3 ) + 3S(0)5(0)5 2 (9, - 9 2 ) + 6 3 (0))J%] , 

where we introduced the notation 

Ft(9^ ...,9 n )= \imF* 2n (\9 n + m + e, . . . ,9, + m + e,0 u . . . , 9 n ) . (12) 

e— +0 ' 

T 2 (9i) is a constant by relativistic invariance. We inserted in Eq. (|T2|) an infinitesimal 
displacent from the pole configurations and will discuss later the result of the limit. For 
the time being we stress that the displacement has to be the same for all rapidities. Any 
different choice would introduce an unjustified asymmetry among the particles. 

Notice that, as a combined effect of the delta functions associated with the discon- 
nected parts, the ^-matrix appears explicitely in the expressions for the /* only through 
the constant factor 5(0). It follows from Eq. (§) that 

S(0) = ±1 . (13) 

Eq. (|10"D shows that S(0) determines whether the statistics of the particles is bosonic or 
fermionic. This means that the expansion of ($>)r in the generic integrable case is formally 
the same than in the free bosonic and fermionic cases, all the difference being encoded in 
the explicit expressions of the form factors. 
The sum of the /* we determined gives 

f(R) = Z { F* + / d9 x e~ ElR [1 + S(0)e~ ElR + e~ 2ElR ] 

2ix J 

+ \ j^y 2 J d9 x d9 2 e- E2R [1 + (e" MRcosh ^ + e - MRcM )S(0)] Ff(9 u 9 2 ) 
+ ^ j^y 3 J d9 x d9 2 d9 3 e-^ R Tt(9 x , 9 2 , 9 3 ) + 0(e' 4MR )} , (14) 
where all the singular delta function terms factorise into the quantity 
Z = 1 + J d9 x e~ ElR 5(0) + - J d9 x d9 2 e"^ [S(0)5 2 (9 1 - 9 2 ) + 5 2 (0)} 

+ i J d9 x d9 2 d9 3 e~ E3R [25(9 X - 9 2 )5(9, - 9 3 )5(9 2 - 9 3 ) + 3S(0)5(0)5 2 (9 1 - 9 2 ) + 5 3 (0)} 
+ 0(e~ AMR ) . (15) 

Since T\ n = 5 Hj o, we have f 1 = Z, so that the singular disconnected parts completely 
cancel out in the ratio (^j) without the need of any specific regularisation. 
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One can continue the computation and work out the next contributions with n > 3 to 
convince himself that the complete result is 



1 



R 



oo 

5 ^ (2^) r 



-RIRcoshBi 



U=l 



•^2* (01> • • • )^n) 



with 



(16) 



(17) 



1 - 1 S(0) e - MRcosh9 ' 

This expression provides the explicit form factor expansion for the one-point function. It 
only remains to check that the specific form factor configuration (O) is not singular due 
to the annihilation poles. For this purpose, take the form factor FQ 2n (\9%, . . . ,6 ) 2ri ) and 
bring the first two particles on the singular configuration 9\ — > 92 + in. This pinching 
produces a pole with the residue given by Eq. (|TT|). We then repeat the operation on the 
second pair of particles by taking #3 — ► 9± + wr, and so on with all the n pairs. This 
procedure produces n poles. The numerator will contain the product 



n 

fc=i 



2n 

n s & 

i=2fc+l 
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2k I 



(18) 



where each string of scattering amplitudes inside the parenthesis is made of factors S(92j + 
^—92k)S{92j—92k)- It follows from (H) and @ that such factors are equal to 1, so that ( fL8|) 
gives n zeros which cancel the poles and leave a finite result. We made these considerations 
on a form factor which differs from fljjP only for the ordering of the particles, namely 
(recall Eq. (|T0D) for an overall phase which does not affect the final conclusion. 

For later convenience, we write down more explicitely the first few terms in the low 
temperature expansion. They read 



R 



I 1 r°° 

= ^ + — [K (r) + S(0)K (2r)] F% + -= / d9 K (2r cosh^) F*(29, 0) 

7T 7T JO 

+ 0(e~ 3r ) } (19) 



with 



MR 



(20) 



For the specific case in which the operator <$>(x) coincides with the trace Q(x) of the 
stress-energy tensor, an alternative and effective way of computing the one-point function 
is provided by the thermodynamic Bethe ansatz approach 0. It is then interesting to 
make a check of the agreement of the two methods in this case. The one-point function 
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of the trace is related to the ground state energy on the cylinder (without the bulk term) 
E{r) as 

(Q)r = (©>fl=oc + 2tt y ±[rE(r)] . (21) 

The TBA allows the determination of E(r) from the knowledge of the S'-matrix in the 
form 

/iif) 
— L(6) coshO, (22) 

L{6) = -S(0) ln(l - S(0)e~ £( - 9) ) , (23) 
where the 'pseudoenergy' e(9) is determined by the integral equation 

e(0) = r cosh 9 - — f d9' <fi(9 - 6')L(6') , (24) 
2n J 

<p(O) = -i±lnS(0). (25) 

The pseudoenergy e(8) tends to r cosh# in the zero temperature limit. Iterating (|4]) once 
with this initial condition and substituting into Eq. (|22| ) produces the result[] 



E(r) 
~M~ 



1 1 2 r°° 
— [ifi(r)+SY0)-Ki(2r)] / d9 K 1 (2r cosh 0)^(20) cosh 0+O(e~ 3r ) . (26) 

7T 2 7T 7 



We can now evaluate (Q)r to this order from (pi] ) using the identity -4- [^^1(2)] 
— -K"o(^)- The final result exactly coincides with (|19[) because it was shown in Ref. 
(Eqs. (3.23), (3.24)) on completely general and independent basis thatQ 



jf = 2vrM 2 , (27) 
Ff{O u 9 2 ) = 8nM 2 ip{9 x - 9 2 ) cosh 2 hzJ!l . (28) 



We also recall here that nearby a fixed point the ground state energy behaves as []10 
E(R) ~ -nC eff /6R, where 

C e// = C - , (29) 

is given in terms of the central charge C and smallest scaling dimension X m i n at the 
conformal point (X min = in a 'unitary' theory). Integration of Eq. (^) then gives 

C ef V f ~ C lff = ^pj drr ^ ~ ^)^) , (30) 



^dn the computation one uses the fact that ip(6) is an even function. 

2 It has been checked for many models that only the e-prescription (|l^) reproduces the general result 
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for the total variation along the whole renormalisation group flow (Cfjj = in a massive 
theory). It is easy to check that inserting in this formula Eq. ([16]) with fl2"TD as the only 
non-zero contribution for n > gives the expected results 1 for the free boson and 1/2 
for the free fermion. One can also check (in the free as in the interacting case) that the 
approximation of the r.h.s. of (|30| ) through the first low temperature terms fll9|) produces 
a poor result, indicating that the convergence to the intermediate and high temperature 
regimes is not very rapid. This fact can be contrasted with the impressively fast con- 
vergence of the form factor expansion (as a function of the distance) for the correlation 
functions in integrable field theory at zero temperature. This different behaviour is hardly 
surprising in view of the different nature of the two expansions. 

Since the distribution (|TT| ) becomes singular as R — > in the bosonic case but not in 
the fermionic case, the high temperature behaviour is very sensitive to the statistics. For 
the quadratic operators in the free cases one has 

vrM 2 /r , S = 1 

[(Q) R ~ (0)i?=oo]/ree = 2M 2 £ S^K^kr) ~ { (31) 



k=l 



-M 2 lnr, S 



For the vertex operators V a = e a{p with scaling dimension X a = —a 2 /4ir, in the free theory 
one has FXn = a 2n /2 n , so that fllfil) gives 

(e atfi )% ee ~ e a2/4r , R^O. (32) 



In the bosonic case S(0) = 1, it follows from the 1/r behaviour of the function (|T7|) that 
the n-particle contribution to behaves as l/r n at high temperature. Up to unlikely 
cancellations of infinities, this means that the integral in the sum rule (^) diverges if (O) r 
receives contributions with n > 1, namely if the theory is interacting. Therefore, one is 
led to the conclusion that in two dimensions the bosonic statistics applies only to the free 
casef] S(6) = 1. An illustration of this fact is given by the sinh-Gordon model, namely the 
theory of a scalar field self-interacting trough the potential fi cosh g<p. Its exact scattering 
amplitude S(6) = (sinh# — i sin 7rl?)/(sinh 9 + isimrB), with B = g 2 /(87c + g 2 ), shows 
that S(0) becomes —1 as soon as g is taken different from zero. One can expect that the 
one-point function (e aip )n will behave in the high temperature limit as a power law with 
an exponent that diverges as g — > in order to reproduce the essential singularity (|32f) . 

We said that the operator $(x) entering (g) has to be local with respect to the excita- 
tions of the zero temperature theory. We now clarify the origin of this requirement. The 
3 See Rcf. |0| for an example of the troubles one encounters when looking for a counterexample. 
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general situation can be illustrated through the example of the Ising field theory defined 
at R = oo by the action 

A = A C FT + r J d 2 xe{x) + h J d 2 xa(x) , (33) 

describing the perturbation of the Ising conformal point by the two relevant operators of 
the theory, the energy e{x) and the spin a(x). 

Consider first the case of vanishing magnetic field, h — 0, in which (^) describes a 
disordered phase (r > 0) and an ordered phase (r < 0) dual to each other and both corre- 
sponding to a free fermionic theory (S = —1). What changes in the scattering description 
of the two phases is the nature of the particles, which are ordinary excitations over the 
unique vacuum at r > 0, and kinks interpolating between the two degenerate vacua of 
the spontaneously broken phase at r < 0. The perturbing operator e is proportional to 
the trace of the stress-energy tensor (G ~ re) and is local with respect to both kind of 
excitations. Its finite temperature one-point function is given by (|3~1~D. The form factors 



of the spin operator are known in both phases [12, 13]. At r > 0, they are non-zero only 
on states with an odd number of particles so that (ct)r vanishes, as it should by spin 
reversal symmetry. In the broken symmetry phase the spin couples to the states with 
an even number of kinks, and one could expect a non-trivial result from Eq. (0). The 
spin, however, is not local with respect to the kinks, and the residue (|TTD on the annihi- 
lation poles gets modified by a constant phase factor (a minus sign in this specific case) 
multiplying the product of scattering amplitudes |I3" |. A similar modification is induced 



in the factor (|T^) which no longer cancels the annihilation poles of the matrix elements 
12), so that the expression (|16D becomes badly divergent for R < oo. The conclusion is 



that (cr)n is not defined at r < 0. Since a finite value of (ct)r would imply the existence 
of spontaneous symmetry breaking on the cylinder, this result illustrates the absence of 
phase transitions in systems with a single infinite dimension. 

The same mechanism applies to the generic case. In two (infinite) dimensions, only 
a discrete symmetry can be broken spontaneously. The elementary excitations in the 
broken phase are kinks interpolating among the discrete vacua. The order parameter is 
non-local with respect to them, and its expectation values on the asympotic states are 
incurably divergent. 

Let us conclude this discussion of the Ising field theory at finite temperature by con- 
sidering the other integrable direction of the action (^), namely the case h ^ 0, r = 0. 
The spectrum of the theory |TJ[ contains now three particles below the lowest two-particle 
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threshold with masses Mi, 

M 2 = (1.6180339887..) Mi, 
M 3 = (1.9890437907..) Mi. 

Hence, the first terms of the low-temperature expansion for the one-point functions are 

(*>*_ = 1 + 1 3 

7T 



_R=oo 



where the ratios^ 



-4 



1 + - Y, Af K (MiR) + 0(e- 2hhR ) 
71 i= i 

<0 1 ^&(0) | a* (iTr) , a< (0)> 



(34) 



(35) 



($>«=co 

are universal and can be extracted from the work of Refs. 1151, IIHI on the form factors of 



the Ising model in a magnetic field. We give their values in Table 1. 



$ 


a 


e 


Af 
Af 
Af 


-8.0999744.. 
-21.206008.. 
-32.045891.. 


-17.893304.. 
-24.946727.. 
-53.679951.. 



Table 1. The universal amplitudes entering the expansion (|34"D for the Ising model in a 
magnetic field. 

Notice that in this case the trace of the stress-energy tensor is proportional to the spin 
operator. Hence, the amplitudes Al are 2irMf divided by (cr)ij =00 . The latter quantity is 
also known from the TBA |l7j. The amplitudes Af, instead, can only be obtained through 
a form factor computation [TJ| in which the 'cluster' property explained in Ref. [18[] plays 
an essential role. 

In conclusion, the form factor approach provides a systematic and explicit low temper- 
ature expansion for the one-point functions in integrable field theory at finite temperature. 
Form factors have already been computed for many integrable theories and it is highly de- 
sirable to have numerical data to compare with the theoretical predictions. It also seems 
natural to expect that the study of the two-point functions should be approachable along 
the same lines with the goal of obtaining a low temperature, large distance expansion. 
This problem is beyond the scope of this note and will be considered elsewhere. 

Acknowledgments: I thank G. Mussardo for interesting discussions. 
4 Here a,i(9) denotes a particle with mass Mi and rapidity 9. 
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